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We investigate the impact of extra vector-like GUT multiplets on the predicted value 
of the strong coupling. We find in particular that Yukawa couplings between such extra 
multiplets and the MSSM Higgs doublets can resolve the familiar two-loop discrepancy 
between the SUSY GUT prediction and the measured value of a 3 . Our analysis highlights 
the advantages of the holomorphic scheme, where the perturbative running of gauge 
couplings is saturated at one loop and further corrections are conveniently described 
in terms of wavefunction renormalization factors. If the gauge couplings as well as the 
extra Yukawas are of 0(1) at the unification scale, the relevant two-loop correction can 
be obtained analytically. However, the effect persists also in the weakly-coupled domain, 
where possible non-perturbative corrections at the GUT scale are under better control. 



1 Introduction 

The consistency of low-energy data with supersymmetric gauge coupling unification p] 
is one of the strongest reasons to expect the discovery of supersymmetry at the LHC. 
Moreover, gauge coupling unification is very well-motivated in heterotic string compact- 
ifications (see [2] for some of the recent developments) as well as in F-theory [3]. Thus, 
if supersymmetry is discovered, the SUSY GUT framework will provide one of the most 
direct ways to access the fundamental high-scale theory (see [I] for a recent phenomeno- 
logical study). 

The SUSY GUT prediction for the strong coupling is, however, not perfect. In fact, 
using two-loop j5 functions and identifying the effective SUSY breaking scale with mz, the 
prediction misses the measured value of otzijnz) by many standard deviations. Possible 
resolutions of this discrepancy include corrections due to an unusual SUSY spectrum 
(for a recent analysis see [5]) or unexpectedly large GUT thresholds. 

In the present paper, we focus on a different possibility: In addition to the usual 
MSSM spectrum, we allow for extra chiral supermultiplets in complete vector-like GUT 
representations. Such states appear in many different contexts and are well-motivated 
theoretically (e.g. as the "messengers" of gauge mediation). While these extra multiplets 
do not affect the one-loop prediction for a 3 , the induced two- loop effect is significant. 
Unfortunately, it further enhances the familiar problems of the MSSM two-loop predic- 
tion [5J. 

However, with extra multiplets naturally come extra Yukawas. We focus on Yukawa 
couplings of the extra multiplets with the MSSM Higgs doublets (in line with the general 
structure of the MSSM and with the natural extension of R-parity to these multiplets). 
As a simple example one may think of a fourth vector-like generation. We also note that, 
following [7HH], a similar class of models has recently been considered in [TUHT2] as a 
possible solution to the little hierarchy problem. Furthermore, extra multiplets may also 
improve the little hierarchy problem in the context of [13]. We take all of this as an 
important extra motivation for our scenario. 

Our results show that the new Yukawa interactions induce a significant shift in the 
strong coupling towards its correct experimental value. Moreover, we find that if the 
extra Yukawas are relatively large at the GUT scale, we end up with an almost perfect 
prediction for a 3 . 

The paper is organized as follows: We start in Sect. [2] with a two-loop analysis of 
MSSM gauge unification in the holomorphic scheme. While gauge couplings run only at 
one loop, the Z factors of chiral multiplets receive corrections at all loop orders. However, 
to achieve two-loop precision for the a 3 prediction, it is sufficient work with one-loop Z 
factors. Their effect on the a 3 prediction comes from the transition to the canonical 
scheme (via the vector and Konishi anomalies), which we perform at the electroweak 
scale. In this approach, the two-loop correction to the MSSM prediction for a 3 arises 
from a sum of terms ~ hi-Z/, where / runs over all flavors, including in particular the 
two Higgs doublets. It becomes apparent that a significant enhancement of the Higgs Z 



factors can provide the desired shift in the a.3 prediction. 

Extra multiplets are introduced in Sect. 13.11 Their Z factors are not important since 
these fields are integrated out above the scale where the transition to the canonical 
scheme is performed. The detrimental effect of extra multiplets on gauge unification 
mentioned earlier arises solely through the increased value of the gauge couplings at 
high energies. The gauge couplings lead to decreasing Z factors as one runs from high 
to low energy scales, and this effect is enhanced in the presence of extra matter. 

In Sect. 13.21 we introduce extra Yukawa couplings. We first consider an analytically 
calculable example with strong gauge coupling at the GUT scale. This can be realized 
introducing extra multiplets in the 10 + 10 and 5 + 5 of SU(5) with masses near the TeV 
scale. Assuming that the extra Yukawas of type 10 10 Hjj and 10 10 Hjj are also strong 
at the GUT scale and neglecting the small effect of the top Yukawa coupling, we solve 
this model analytically. The resulting shift in the a 3 prediction, which is due to enhanced 
Higgs Z factors, leads to nearly perfect agreement with the experimental value. 

In Sect. 13.31 we extend our analysis to more general scenarios, keeping in particular 
the GUT-scale gauge coupling in the perturbative domain. We demonstrate that our 
promising initial results retain their validity in this setting as long as the extra Yukawa 
couplings at the unification scale are sufficiently large. 

Scenarios with a larger number of extra Yukawa couplings are investigated in 
Sect. 13.41 We focus on two particular types of extension. First we analyze models in 
which further 10 + 10 pairs with further extra Yukawas are introduced. We find that in 
this case the total effect on the Higgs Z factor can not be increased significantly. The 
second type of models we consider possess the same matter content as our minimal model 
from Sect. 13.31 (i.e. a 10 + 10 plus a 5 + 5 pair). However, we now also allow for couplings 
of type 5 10 Hp and 5 10 Hy. We find that in this case the two-loop prediction for the 
inverse coupling 2ir/a 3 is increased even further. In particular, we are able to reproduce 
the 'optimal' results from Sect. 13.31 under milder assumptions (i.e. for lower couplings at 
the GUT scale). 

Low- and high-scale threshold corrections and, in particular, the critical issue of 
strong gauge couplings at the GUT scale [T4Tl2~T] are discussed in Sect. HI We emphasize 
that, due to the absence of higher-order perturbative corrections to the holomorphic 
couplings, precision unification is not compromised when allowing for relatively large 
values of the GUT coupling. However, once one reaches the actual strong-coupling regime, 
non-perturbative corrections can arise and affect the a 3 prediction. Thus, it appears to 
be safe to stay in a domain where terms that are exponentially suppressed by the inverse 
gauge coupling are negligible. This is sufficient for our purposes. 



2 Conventional unification from a 
holomorphic perspective 

2.1 Basic formulae with holomorphic thresholds 

We find it convenient to work in a holomorphic scheme, where the perturbative running 
of gauge couplings is saturated at one loop [22H24] . Focussing for simplicity on SU(5) 
models with adjoint breaking, we have 

2vr 2tt , , M , f3l , M , XYl M - M 

+ h\n — + 6 3) ln^ + 6f' y ln^— + 6f In -=- , 1 
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where a>h,i(n), with z = 1,2,3, are the holomorphic U(l), SU(2) and SU(3) gauge cou- 
plings of the MSSM at some low scale ft and a^(M) is the unified holomorphic coupling 
at some high scale M. This equation is exact to all orders of perturbation theory provided 
that fj, lies above the supersymmetry breaking scale. While the first logarithm combines 
the effects of all fields that remain light below the GUT scale, the last three logarithms 
are associated with heavy fields. The various fej's are the appropriate one-loop /3-function 
coefficients. They are labelled by '(3)' for Higgs-triplet, l X, V for the massive vector 
multiplets of X, Y gauge bosons, and $ for those components of the GUT-breaking field 
$ that are not 'eaten' by the X,Y bosons. We emphasize that mj, m^ and m| are 
holomorphic rather than physical mass scales. In other words, these are the mass param- 
eters of the holomorphic Wilsonian action, in which kinetic terms are not canonically 
normalized |2] 



At the moment, we may think of a situation where va\ ~ m^ y ~ m| ~ -^gut an d 
fi <C Mq UT <C M. Furthermore, we identify the SUSY-breaking scale with mz for now, 
postponing a brief discussion of low-scale threshold effects to Sect. HI 

Thus, we can set /x = mz and translate holomorphic to canonical gauge couplings 
using the well-known anomaly relation [22Tl24) 1 l : 

2tt 2tt 
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Here c^ = gf /(47r) is the canonical gauge coupling, T(Gi) = C 2 (Crj) is the Dynkin index 
or quadratic Casimir of the adjoint of the gauge group G^ and T(RJ) is the Dynkin 
index of the representation -R/of the flavor / (with respect to Gi). This gives rise to the 
'holomorphic master formulao 
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(3) 



*At our level of accuracy, we can ignore corrections associated with the transition to the scheme- 
dependent (e.g. DRED) physical gauge coupling [25] . 

2 i.e. the holomorphic version of what is called the 'master formula' in [l6l[23] 



Obviously, the choice of M in this formula is irrelevant for the prediction of a.3 from a.\ 
and a 2 - In particular, we can set M = m XY (irrespective of the actual UV completion 
scale) and write 

-^4 = (^ S + bi ln — - T ^ X <9 2 i{mz)) - E T (^/) foZfimz) + Af, GUT , 
ai{m z ) a h {m XY ) m z ^ 

(4) 
with (holomorphic) GUT threshold corrections 



AjW^In^+tfln^- (5) 



The above threshold corrections will be small if the two relevant mass ratios are 0(1). 
This will be the case if the superpotential contains no parametrically small couplings. 
Indeed, staying strictly within the holomorphic scheme, not only m\ and m^, but also 
m x Y are determined purely by superpotential terms (in contrast to the physical masses 
of the X, Y multiplets, which include a prefactor g coming form the gauge kinetic term). 



2.2 Compatibility with the conventional master formula 

In order to establish the equivalence between Eq. ([3]) and the more conventional master 
formula of [23], we need to replace holomorphic by physical mass parameters: 

mP XY = gm h XY Z l J 2 , m%Z^ = m% , m p 3 Z 3 = m% , (6) 

where g, Z§ and Z 3 are the (GUT-scale) gauge coupling and Z factorsj and the super- 
script 'p' stands for 'physical'. Using these relations together with Eq. ([2j), we rewrite 
Eq. (J3J) in terms of canonical gauge coupling and physical masses. Making also use of 
the identity bf' Y - 2T(Gi) = -2T(SU(5)), which follows from bf' Y = -2T(R?' Y ) and 
T{Rf' Y ) + T(Gi) = T(SU(5)), and choosing M = m XY , we find 

27r 2tt m XY ai(m z 



+ bi ln H^L _ T(Gi ) , n *W _ £ T{R{) ,„ ^M 



ai{m z ) a{m XY ) ' m z ' " ai(m p XY ) ^f^' ' ' '" Z f (m p * 

with (physical) GUT threshold corrections 



^p 

X,Y> ~T ^f\" b X,Yi 

(7) 



A P _ b m ln m X,Y b * ln m X,Y 

^i,GUT - °i m r> + °i m v 

"t>3 '"'$ 



This is in agreement with |6][T6][23] . It simply represents a slightly different parame- 
terization of our ignorance of the high-scale input: In the 'holomorphic master formula' 



3 Since they appear only in logarithms, it does not matter whether we use a holomorphic or a 
canonical gauge coupling and at which precise mass scale we evaluate all these quantities. 



of the previous subsection, all the relevant non-holomorphic (and in this sense 'unpro- 
tected') input data were the high-scale boundary values of the Higgs Z factors. (Note 
that high-scale Z factors of complete SU(5) multiplets, such as MSSM matter, do not 
affect the a 3 prediction.) Now, by contrast, the very same ignorance is hidden in the 
value of the physical Higgs triplet mass. 



2.3 Predicting a% 



It will be convenient to rewrite Eqs. (j4]) and (JSJ) as: 
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collecting all two-loop effects (the logs of gauge couplings and Z-factors as well as GUT 

thresholds) in the correction terms Aj. Multiplying the first of these three equations by 

ip 2 — b 3 )/(bi — b 2 ) = 5/7, the second by (b 3 — b\)/{b\ — b 2 ) = —12/7, and adding them to 

the third equation, one finds 

->„. 1 o o^n r k 1 o 
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one— loop prediction 




two— loop corrections 



Here we have suppressed the mass scale, mz, for brevity. Setting the holomorphic 
GUT thresholds to zero, the two-loop correction, i.e. the second bracket in Eq. (fTU|) . 
explicitly reads 
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Given that we are only aiming at two-loop accuracy, it is consistent to evaluate all 
quantities entering the above expression with one-loop precision. In particular, we can 
use one-loop Z-factors and the experimental values for the gauge couplings (given by 
2ir/ai = 370.7, 2ir ja 2 = 185.8, 27r/«3 = 53.2). The Z factors can be written as 



J L,Hu,H D 

Zq 
Ze 
Zu 
Z D 



J SU(2) 



x Z, 



Zsu{z) x Z 

Z f/(1) 



U(l) 
SU(2) 

16/9 



X Z 



1/9 

U(l) 
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Introducing the shorthand notation Zjjm 
tive gauge group contributions read 



SU(3) X Z u(1) . 

~ Z\i Z S u(2) = Z 2 , 



'SU(3) 



Z 3 , the respec- 
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Here Cp,i are the quadratic Casimir operators of the fundamental representations of 
U(l), SU(2) and SU(3), whereas bi are the respective one-loop /3-function coefficients (cf. 
App. A). Furthermore, ckgut is the one-loop GUT coupling, 27r/a GU T = 153. Equation 
( Ti"3"j) is valid under the assumption that all Z-factors are unity at Mqut- 

Combining all contributions, we find: 

^ = 53.7- 4.08 -0.64-5.51 + 0.21 + 3.22 + 1.83- 0.42 

1— loop vector anomaly Z^ Zq Z^ Zjj Zp Zfj ,Zfj 

= 53.7 -5.4 = 48.3, (14) 

i.e., the two- loop corrections shift the one- loop prediction away from the experimental 
value of 53.2. As already advertised in the Introduction, this set of numbers suggests a 
simple way to cure the problem: It will be sufficient to introduce a significant enhance- 
ment of the Higgs Z factors at low energies. This will be realized in the following using 
extra multiplets with extra Yukawa couplings. 

For completeness, we record the result which is obtained if the canonical (rather 
than the holomorphic) GUT threshold corrections are assumed to vanish. At the technical 
level, this corresponds simply to replacing the 'vector anomaly' contribution —4.08 above 
with 

24 k ^N _ 3h 03N = _ 06? _ 317 = _ 3 ^ (15) 

7 CtGUT «GUT 

The resulting prediction improves insignificantly. 

The smallness of this change is due to the weak dependence of the a% prediction 
on the value of «gut in Eq. ( TT5l) . This is the result of the approximate cancellation 
24/7 — 3 = 3/7. It is equivalent to the statement that the splitting between itlx,y and 
m<j> affects the a^ prediction only very weakly, which reflects the similarity of the ratios 
of the bf and the ratios of the corresponding MSSM coefficients &,. 

Finally, we comment on the effect of MSSM Yukawa couplings which we have so far 
neglected. At moderate tan/3, the top Yukawa coupling gives the dominant correction 
(cf. App. B), 

Atop (^) =0 ' 32 ' (16) 

which is however far too small to cure the 'two-loop' problem. 



3 Extra multiplets 

3.1 Extra multiplets without Yukawa couplings 

To preserve one-loop gauge coupling unification in the most straightforward way u, we 
restrict our attention to complete SU(5) multiplets. More specifically, we focus on models 



l See [26 for an alternative point of view 



with n 5 pairs of 5 + 5 and n w pairs of 10 + 10. At one loop, this leads to a modification, 
bi — >■ b\ = hi + n, of the MSSM /3 function coefficients, where n = n 5 + 3n w is known as 
the 'messenger index'. This, of course, does not affect the one-loop a% prediction. 

The two- loop correction, given explicitly in Eq. ([lip , changes only because of the 
modified Z factors of the MSSM matter fields. We emphasize that Eq. (1111) does not 
need to be supplemented with Z factors of the extra multiplets since these are assumed 
to decouple above mz- 

The modified Z factors are obtained from Eqs. (fT2l) . as before, but now with the 
gauge group contributions (cf. Eq. (TT3|) ) 



Zdmz) 
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Here o-gut is the one-loop GUT-coupling in the presence of extra multiplets. Note that 
the one-loop GUT scale, Mqut = 2 ■ 10 16 GeV, is unaffected by the presence of the 
additional matter. For simplicity, we have at his stage neglected the necessary hierarchy 
between mz and the scale at which the extra multiplets decouple. 

The two-loop corrections related to the different matter field Z-factors can be found 
in App. C. Here we only list the overall two-loop prediction for 27r/a 3 for different values 
of the n parameter: 

Table 1: Two- loop prediction for 211/0,3 



n 


Prediction for 2ir ja^ 


1 


48.08 


2 


47.91 


3 


47.61 


4 


46.81 


4.45 


45.81 



The value n = 4.45 formally corresponds to Squt = 1. Of course, this has to be 
interpreted as n = 5 together with an appropriately raised decoupling scale. We see that 
with increasing n the two-loop prediction for a 3 becomes systematically worse. 

Of course, we could easily repeat the analysis using the assumption of vanishing 
canonical GUT-thresholds. In this case, the increased value of the GUT-coupling enters 
the result also via the analogue of Eq. ( TT5"|) . with ckqut replaced by «gut- For n = 4.45, 
this contribution changes from —3.84 to —5.16, giving 2ir /0J3 = 44.73. 

3.2 Extra multiplets with Yukawa couplings - an example with 
strong GUT coupling 

In the following we extend our analysis of Sect. 13 .11 by allowing renormalizable couplings 
between the SU(5) matter and the MSSM particles. To this end let us assume that there 



8 



is at least one pair of 10 + 10 vector-like matter. In analogy to the MSSM, we denote 
its field content by 10 = (Q e , U e , E e ) and 10 = (Q e , U e , E e ), with an index 'e' for 'extra 
multiplet'. We can now introduce extra Yukawa couplings of the form: 

WD nQ e U e H v + K~Q e U e H D , (18) 

These new interactions modify the prediction for a% solely through their effect on the 
Higgs wavefunction renormalization factor Zh (where Zh stands for either Zh v or Zh d ). 
We note that the couplings k, k can be extended to full 4x4 Yukawa matrices allowing 
for mixing between the three MSSM generations and the additional matter. However, 
these mixings have to be small due to FCNC constraints [11] (see also [57|BB]) and we 
neglect them. 

In order to simplify the analysis in this section let us once again neglect any contri- 
butions arising from the MSSM Yukawa sector (based on our results from Sect. 12.31 and 
App. B we expect those contributions to be small). This allows us to treat k and 75 on an 
equal footing. In particular we can set k = K and effectively deal with a single Yukawa 
coupling (say k) and a single Higgs (calling this field H). 

The one-loop RGE for Zh (with t = ha.fi) reads 

d\nZ H 3 3 

2tt — — — = -3a K + — o-i + -a 2 , (19) 

at 10 2 

where we have defined a K = k^/(47t), with k 2 = k, 2 /{ZhZq^u,,) the canonical Yukawa 
coupling. It is clear that a large k will drive Z H to larger values at the electroweak scale, 
improving the a^ prediction. 

Equation ( II 9p entails the factorization property 

Z H = Z% x Z Y H , (20) 

gauge part Yukawa part 

where Z H and Z H represent the contributions from the gauge couplings ai, a 2 and the 
Yukawa coupling a K . As before, the gauge part Z§ is determined by Eqs. (TT2"j) and ( TTTj) . 
We thus focus on the one- loop running of Z^: 

The corresponding one-loop RGE for a K reads 

dhia K 16 13 

2vr — -— = 6a K -a 3 - 3a 2 - — «i • (22) 

at 3 15 

In the following we will analyze a model in which the value of the low-scale Zj^-factor 
can be obtained in a completely analytical manner. To this end let us assume that both 
the extra Yukawa couplings as well as the gauge couplings begin their evolution at the 
strong-coupling point at the high scale: 

a GUT ~ a K (M GUT ) ~ 1 . (23) 



Formally this corresponds to n — 4.45. In this case the relations 

U b' 

a 2 = 77^3 «i = -f a 3 (24) 

will be approximately valid at all energies significantly below the GUT scalejj Equa- 
tion fl22|) now takes the form 

dlna K „ {16 3b' 3 13b' 3 



2 *^ir = 6a *- a °{-3 + ii + mJ- (25 > 



Using the one-loop RGEs for the gauge couplings 



^=-»- ,26) 

this can be further rewritten as 

d\n(a K /a 3 ) „ /16 3b\ 1364 „\ 

Eq. ( )27l) has an infrared-stable fixed point of the Pendleton- Ross type [29] given by 

a K = 1.37a 3 for n = 4.45 . (28) 

As a result of fast initial evolution, the ratios of Yukawa and gauge couplings quickly 
reach the fixed-point regime, which is then maintained all the way down to the weak 
scale. From Eq. (12 lj) we now have 

dlnZYr _ a 3 4.11 dlna^ 1 , , 

— - = -3 • 1.37 — = — -2- 29 

dt 2n b' 3 dt v ; 

and 

In Z Y H (m z ) = 2.83 \n(a GVT /a 3 (m z )) = 6.05 . (30) 

The resulting correction to the 0:3 prediction is 

&J— )=l In Z%(m z ) = 7. 78, (31) 

V a 3/ 7 

which is just sufficient to compensate the negative two-loop effects in the last line of 
Table [TJ Of course, at this stage our promising results should be taken with caution. 
Specifically, when talking about a strongly-coupled unified theory, one faces the danger of 
potentially large and incalculable corrections at the GUT scale which could, in principle, 
render the entire two-loop analysis obsolete. We postpone the discussion of these issues 
to Sect. SJ We also note that the influence of extra Yukawas above the GUT scale on the 
unified coupling has recently been discussed in [30] . 



One writes a^" 1 ^) = 0(1) + 6- In (^^j and neglects the 0(1) term. 

10 



3.3 Effect of extra Yukawas in models with perturbative gauge 
couplings 

In this section we extend our previous results to a more general setting by lifting some 
of the simplifying ad hoc assumptions that were made so far. This means in particular 
that we introduce an explicit decoupling scale M for the extra matter fields. Also, from 
now on the messenger index n is allowed to attain only integer values (however, we only 
allow for n > 3 since we assume at least one 10 + 10 pair). The superpotential in the 
extended Yukawa sector is still specified by Eq. (1TB]) . In order to improve the accuracy 
of our predictions we will also take into account the contributions from the MSSM top 
Yukawa coupling. Since this coupling enters the one-loop RGEs of a K and a-^ in a different 
manner, this step explicitly breaks the symmetry of our model with respect to k and K. 
In particular, from now on we will distinguish explicitly between these two couplings. 

The most important constraint which we impose on the models in this section is 
perturbativity of the gauge couplings. Among other things this implies that the analyt- 
ical approach we developed in the previous section is no longer applicable. Instead we 
will employ an alternative technique, which will allow us to handle our models semi- 
analytically. 

To this end we will have to deal with the two different regimes of the theory (the 
high-energy regime above M and the low-energy regime below M) separately. Let us 
focus on the high-energy theory first. As was mentioned before, any MSSM field obeys 
a factorization property analogous to Eq. ( 120]) . This means that the gauge and Yukawa 
contributions to the matter field Z-factors decouple and can be analyzed independently. 
The gauge parts of the Z-factors are obtained from Eq. ( TT71) by simply replacing rriz — > 
M. In order to get the Yukawa parts we will need the one-loop RGEs for the Z Y (above 
M): 

d In Zjj d In Zjj 

2n — - = — 3a t — 3a K , 2tt — — = — 3«k , 

dt ' dt 

d\nZ Y dlnZX / N 

2vr — ^ = -2a t , 2vr -^ = -a t . 32 

dt ' dt K J 

Since we regard only the k, k and yt Yukawas as non-vanishing, all other Z y -factors are 
irrelevant for our analysis. In the following we will also need the one-loop equations for 
the three aforementioned Yukawa couplings: 

dlna t 13 16 

27r — -— = <6a t + 3 a K - — a\ - 3a 2 — a 3 

dt 15 3 

dhia K 13 16 

2vr — — — = 6a K + 3a t - — a x - 3a 2 — a 3 (33) 

dt 15 3 

d\na K 13 16 

27r — — — = 6a^ - -— «i - 3a 2 — a 3 

dt 15 3 

Note that these identities generalize Eq. ( )22l) to the case of non- vanishing at- Following [6] 
we can combine Eqs. (133]) and (132]) to get a closed analytic expression for the Yukawa 
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parts of the relevant Z-factors: 
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The theory below the scale M is the MSSM. Therefore we can apply our formulas 
from Sect. 12.31 and App. B. The only modification is that we now integrate from m z to 
M rather than to Mqut- As an illustrative example we record the result for the Yukawa 
part of the Zh u factor: 



Z H v ( m z) 

ZUM) 



a 2 



«i(M) \ 30i, i 

Oii(m z )J \a 2 (m z )J 



(M) \ 



3 
2 6., 



q 3 (M) \ 

a-i{m z ) ) 



8 
3 63 



ott{M) \ 
oit{m z )) 



(35) 



The calculation of the other Z-factors proceeds in a similar manner. The low-scale values 
Z(m z ) are then obtained by multiplying the expressions for Z(M) and Z(m z )/Z(M). 



The brackets in Eqs. (I34|) and (135]) involving gauge couplings can be evaluated ana- 
lytically (by using the respective one- loop values). To calculate the Yukawa brackets we 
have solved the one-loop RGEs in Eq. (1331) numerically by evolving them from the GUT 
down to the electroweak scale and using a ^-function approximation at the decoupling 
scale M. The resulting low-scale values for the three Yukawa couplings were then substi- 
tuted in Eqs. (l34"j) and (1351) . In Tables [2] and [3] we have listed the two- loop corrections to 
27r/a!3 for different values of the two parameters M and n. We have also tested the sen- 
sitivity of our results against variations of the initial values a K (M GU T) and o^Mgut) 
(we remark that the values « k (Mgut) = «k(A^gut) = 0.228,0.457,0.913 correspond to 
1/6, 1/3, 2/3 times the fixed point value 1.37 of the extra Yukawa couplings). In each case 
the initial value for the top Yukawa coupling at Mqut has been adjusted to reproduce 
the correct low-energy parameter y t (m z ) = 0.99 (see also App. B). 

A quick glance at Tables [2] and [3] reveals that models with n = 4 are favoured over 
their n = 5 counterparts. The former are also of particular interest because of the low 
value of M [TUHl"2"| l3"l"] . We have intentionally omitted the n = 3 case because the n = 3 
models are unable to generate a sufficiently large GUT coupling (say «qut > 0.2). 

It is important to note that the low-energy couplings a K {m z ) and a- K {m z ) are virtu- 
ally insensitive to their input values at the GUT scale. This observation indicates a very 
straightforward way of increasing the two-loop prediction for 2^/0:3 - namely by taking 
the input parameters a Kj GUT and o^gut as large as possible. Note also that (in con- 
tradistinction to gauge couplings) the Yukawas do not exhibit any flavor enhancement 
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Table 2: Numerical results n = 4 (one pair of extra Yukawa couplings) 



M(GeV) 


acuT 


Ok,GUT — Ok.GUT 


a K (M) 


<*k(M) 


y t (M GVT ) 


2n/a 3 


500 


0.227 


0.228 


0.102 


0.134 


0.68 


51.10 






0.457 


0.102 


0.134 


0.84 


51.55 






0.913 


0.103 


0.134 


1.01 


51.99 






2.000 


0.103 


0.134 


1.26 


52.50 






4.000 


0.103 


0.134 


1.52 


52.95 






6.000 


0.103 


0.134 


1.70 


53.21 


300 


0.245 


0.457 


0.105 


0.137 


0.78 


51.67 






2.000 


0.105 


0.137 


1.18 


52.62 






4.000 


0.105 


0.137 


1.42 


53.07 




0.245 


6.000 


0.105 


0.137 


1.58 


53.33 


Table 3: Numerical results n = 5 (one pair of extra Yukawa couplings) 


M(GeV) 


acuT 


«k,GUT — Ok.GUT 


<*k(M) 


a«(M) 


y t (M GVT ) 


2ir/a 3 


25 • 10 4 


0.229 


0.228 


0.097 


0.121 


0.70 


50.34 






0.457 


0.097 


0.122 


0.86 


50.79 






0.913 


0.098 


0.122 


1.04 


51.24 






2.000 


0.098 


0.122 


1.29 


51.74 


O 4 




4.000 


0.098 


0.122 


1.56 


52.19 


o 4 




6.000 


0.098 


0.122 


1.73 


52.45 



(cf. Sect. Hj). Therefore the actual expansion parameters are a K /47r and a^/Air, which 
means that the strong coupling condition reads a K /47r ~ a-^/An ~ 1. Following this line 
of thought we have considered models whose input values a K; GUT and oj^gut are as high 
as 6.0. 

We once again emphasize that the gauge couplings in this section are only allowed to 
attain perturbative values (in contradistinction to their Yukawa counterparts). As will 
be argued in Sect. HJthe lowering of the unified coupling from «gut ~ C(l) down to 
~ 0.2 can potentially have a dramatic impact on the magnitude and calculability of the 
high-scale threshold corrections around the GUT scale. 



3.4 Models with further Yukawa couplings 

In the following we will investigate the principal effect of introducing further non- 
standard Yukawa couplings. We focus on two extensions: First we consider adding new 
10 + 10 pairs which couple to the observable sector through an interaction analogous 
to Eq. (fl8|) . In particular we increase the messenger index to n > 6. Second, we ex- 
plore a different possibility by restricting ourselves to the n = 4 case but allowing for 
renormalizable interactions between the 10 + 10 and the 5 + 5 fields. 

Following this line of thought let us now increase the number of additional vector-like 
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multiplets to n = 6 by introducing two 10 + 10 pairs. We postulate a superpotential of 
the form: 



W D KQ e U e Hu + nQ e U e H D + k! Q' e U' e Hu + n'Q e U e H D 



(36) 



where the primed fields denote the matter content of the new 10 + 10 pair. In the 
following we will treat the two Yukawa pairs (k,k) and (k',~R?) on an equal footing, i.e. 
we assume that k — k' and re = ~k' . Note that by going from a superpotential of the form 
( 1T81) to superpotential of the form ( 136]) we only change the running of the Yukawa part 
of the Z-factors above the decoupling scale. The relevant one-loop RGEs are 



2ty 
2ty 



dhiZl 

H U 

alt 
dlnZ% 3 

dt 
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-2a t , 
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2n 
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dt 
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dt 



-6«7c- 



-a t 



(37) 



for the Z Y and 
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2n 
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dt 
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9a K + 3a t - — a 1 
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Q 13 

15 
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, 16 
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o 

16 



(38) 



for the Yukawas. 

Combining Eqs. ( 1371) and ( 1381) we arrive at the analogue of Eq. (|34 

zU M ) = 

Zl(M) = 
Z Y Hd {M) = 
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The results from our numerical analysis are listed in Table HI The prediction for a 3 
improves only slightly in comparison to the 'optimal' n = 4 models with a single pair of 
extra Yukawas. 

The next type of models we consider contain one pair of 10 + 10 and one pair of 
5 + 5 extra multiplets. Using the standard decomposition 5 = (D e , L e ) and 5 = (D e , L e ) 
we introduce a superpotential of the form: 



W D KQ e U e H v + KQ e U e H D + XQ e D e H D + XQ e D e H, 



u 



(40) 
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Table 4: Numerical results n=6 (two pairs of extra Yukawa couplings) 



M(GeV) 


«GUT 


«k,GUT — Ok.GUT 


a K (M) 


<*k(M) 


y t (M GVT ) 


2n/a 3 


17 10 6 


0.227 


0.228 


0.064 


0.080 


1.26 


51.08 






0.457 


0.065 


0.080 


1.68 


51.66 






0.913 


0.065 


0.081 


2.24 


52.25 






2.000 


0.065 


0.081 


3.08 


52.92 






4.000 


0.065 


0.081 


4.06 


53.52 






6.000 


0.065 


0.081 


4.77 


53.86 



Note that this is a direct extension of the n = 4 models from Sect. 13.31 - we have simply 
added two new interactions to our superpotential. The calculation of the two-loop as- 
correction proceeds exactly as before. Here we only list the relevant one-loop RGEs above 
the decoupling scale M: 



2ty 
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dhxZl 
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dt 
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— 3«t — 3a K — 3aj , 
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2ty 
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(41) 



for the matter field Z y -factors and 
din at 
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15 3 
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(42) 



—2 



for the Yukawa couplings. As before we have defined a\ = A 2 /(47r) and a j = A /(47r) . 
The modification of Eq.( |34"l) reads 
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(43) 
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We have presented the results from our numerical analysis in Table [51 It is impor- 
tant to note that, already for moderate input values of the extra Yukawas at the GUT 
scale, we reach the region of the experimentally measured a 3 . Hence in models of this 
type it is no longer necessary to invoke excessively large Yukawa couplings at the high- 
scale. For this reason we have restricted ourselves to a region of parameter space where 

«k,GUT, Ok,GUT, OA,GUT, «A,GUT < 1- 

Table 5: Numerical results n=4 (two pairs of extra Yukawa couplings) 



M(GeV) 


CtGUT 


«k,GUT — «k,GUT — OA.GUT — a A,GUT 


yt(M GUT ) 


2n/a 3 


1000 


0.206 


0.228 


2.25 


52.39 






0.457 


3.38 


52.99 






0.913 


5.20 


53.58 



4 Threshold corrections and higher-order effects 



In this section we discuss issues related to low- and high-energy thresholds and other 
higher-order corrections. The effect of the superpartner spectrum on the value of the 
strong coupling has been studied in detail in [32TI34"] . The analysis reveals that the low- 
energy thresholds potentially shift the predicted value of a 3 by a significant amount. For 
instance, it is well-known that certain SUSY spectra with light gluinos can compensate 
for the detrimental two- loop effect discussed in Sect. 12.31 and therefore bring the a 3 - 
prediction in line with the experimental value (see e.g. [US])- However, gluinos tend to 
be heavy in the simplest mediation scenarios and it generally requires a compensation 
of several mediation effects to make them light. A detailed study of concrete models 
realizing this possibility has recently appeared in [5]. 

Another potentially important contribution arises from heavy particle thresholds. 
For example, from Eqs. (J5J) or flS} it is clear that in order to shift the prediction for 
2ti /a.3 by several units the logarithms of the mass ratios itixy/ 171 ^ or m x y/ m 3 have to 
be several units themselves. In other words, the Higgs triplets have to be ~ 10 2 lighter 
than Mgut and proton decay has to be avoided through some version of the missing 
partner mechanism (see [36] for references). Of course, thresholds with larger numerical 
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prefactors (and hence smaller required mass ratios) can arise in models with large GUT- 
scale representations (see, e.g., [32|l37H39j). A similar enhancement can come from large 
multiplicities of heavy states (this has in particular been argued in the context of certain 
string-motivated models |in])Jj All of this clearly makes 'GUT-scale thresholds' a viable 
explanation of the precise value of a^. 

The models considered in this paper offer an alternative solution to the two-loop 0:3- 
discrepancy. This solution differs conceptually from the aforementioned approaches as it 
does not rely on any type of threshold effects. It realizes a lower 0:3 value at the expense 
of relatively large extra Yukawa couplings. The latter do not require any additional 
SU(5)-breaking effect (beyond the doublet-triplet splitting, which is anyway present in 
the MSSM). Nevertheless, precision at the high scale remains a critical issue and the rest 
of this section is devoted to its analysis. 

The extra multiplets supporting the extra Yukawas raise the value of the GUT-scale 
gauge coupling. It is then tempting to consider the extreme case of such scenarios: Grand 
unification at the strong-coupling point. Even more conservatively, one might want to 
drop the GUT-assumption altogether and to demand only that all three SM gauge factors 
become strongly coupled at the same energy scale [H] (for early related work see [H]). 
The resulting high-scale error for the 03 prediction can be estimated using the familiar 
one-loop formula 

+ 2h ln(M GUT /m z ) . (44) 



Oi(mz) a>i(M GVT 

Naively, 'strong coupling' means that the loop-expansion parameter g 2 /(16ir 2 ) is of order 
one. This would imply that 47r/aj(MGUT) = 1 ± C(l) m Eq. f|44l) . The resulting error of 
the «3 prediction is rather small, around 1%. However, because of the large number of 
flavors, the actual expansion parameter is in fact closer to a rather than a/(47r)Jj Being 
at strong coupling then means that l/aj(M GUT ) = 1 ± 0(1), which corresponds to a 
~ 10% error of the a% prediction. While this easily brings the 2-loop MSSM prediction 
for «3 in line with the data, it also makes any discussion of 2-loop effects obsolete: The 
error is simply too large. 

For the purpose of the present paper, we adopt a different point of view: We assume 
that a model with true, calculable unification exists in principle and that the coupling 
strength is controlled by some high-scale holomorphic parameter (e.g. a string theory 
modulus). When talking about strongly-coupled unification, we assume that this model 
is realized in a region of its parameter space where «gut — 0{1). 



6 For explicit orbifold constructions where the relevant corrections decouple from the string scale and 
their size could be checked straightforwardly see e.g. |32]. 

7 There are at least two ways to see this: First, we focus on the contribution of all 'flavors' to the 
one-loop /3-function of the QCD coupling at high scales. The corresponding /3-function coefficient is 
2&| avor = 2(n + 6) w 20 (n — 4 or 5). This more than compensates for the suppression by 4n, leaving 
us with a number close to a as the actual expansion parameter. Alternatively, we consider the one-loop 
contribution of the gluons/gluinos, 2fo§ olor = 18. Once again, this is more than sufficient to cancel the 
factor of 4-7T. 
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To be more specific, consider the string-theoretic (heterotic) formula [43] 

fi(S,T) = k t S + Ai(T) (45) 

for the gauge- kinetic functions of the three SM gauge groups. For the conventional embed- 
ding of Gsm in SU(5) and SU(5) in E 8 , the ki are unity, corresponding to tree-level uni- 
fication. The crucial point is that the modulus governing the tree-level coupling strength 
(in this case the dilaton superfield S) does not appear in the loop corrections. They de- 
pend on a set of different moduli which we collectively denote by T. The reason for this 
is basically the same as in field-theoretic arguments for 1-loop running: holomorphicity 
and shift symmetry of Im(S') [22j[23]. Now, moving the modulus S from its perturbative 
value to the region where the expansion parameter is 0(1), we see that the high-scale 
non-universal correction Aj(T) is not enhancedo The only potential danger comes from 
non-perturbative extra terms ~ Cjexp(— diS) which can appear on the r.h. side of Eq. 
f )45|) . To keep such terms under control, we only have to assume that exp(— ReS) <Cl- 
a much weaker requirement than 1/ReS <C 1. 

A further important issue is the error which builds up along the RG trajectory from 
m z to M GU t as a result of using the two-loop instead of the full /3-function for the gauge 
couplings. According to the previously discussed master formulae (cf. Eqs. (J3J) and (J7|)) 
this error is identical to the error of the In Z terms. At one loop (and focussing only on 
the gauge sector for simplicity) we have d In Z/dt ~ a , which upon integration gives 
In Z ~ Int. It is clear that this leading order effect receives contributions from the entire 
integration range. We will now argue that higher-order corrections are UV dominated. 
To this end we recall that at two loops the RGE for a generic Z-factor has the form: 

d\\\Z 9 /d . 

— a + a 1 , (46) 

where a is the canonical (or physical) gauge coupling. From the anomaly relation Eq. ((2|) 
we obtain schematically 

a ~ a h + a 2 h \na h + a 2 h In Z ~ - + — , (47) 

t t z 

where we have only displayed the leading corrections. Note that we have used the (per- 
turbatively exact) one-loop holomorphic gauge coupling ah ~ 1/t as well as the fact that 
(at one loop) In Z ~ In t . Substituting the above relation in the two-loop RGE f )46|) gives 

, „ f dt f hit , , .„ s 

\nZ~ J - + J — dt, (48) 

where terms ~ 1/t 3 or higher were neglected. The first term on the r.h. side gives 
the previously discussed leading order hit effect. The second integral is UV dominated, 



8 The dual situation, where the leading-order gauge coupling is governed by the GUT-brane volume 
T and corrections (related to a higher-dimension operator) depend on the dilaton S, arises in F-theory 
GUTs. These corrections tend to aggravate the two- loop discrepancy for a^ [31], potentially making 
'our' Yukawa effect the more interesting (see also [45]). 
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i.e., this subleading effect can indeed be neglected at our level of accuracy. To be more 
precise, the corresponding correction is not enhanced by the parametrically large quantity 
t ~ ln(MQ\j^/mz) or a log thereof. Hence, it corresponds to an 0(1) correction to \nZ. 
This is equivalent to a multiplicative 0(1) uncertainty of the high-scale Z factor, which 
we anyway have to accept in the absence of an explicit GUT model. The same argument 
goes through for contributions of even higher order. 

To summarize, our strongly-coupled unification scenario is defined as follows: The 
unified gauge coupling is taken to be relatively large, ogut ~ ®(X)i while non- 
perturbative corrections are still under control. This may be the case because there is at 
least a (small) hierarchy of the type exp(— 47r/acuT) < 1 or because the coefficients of 
such non-perturbative terms happen to be small. In such a setting, our 2-loop analysis 
of the strongly coupled model from Sect. 13.21 is meaningful and necessary. 

We also emphasize that the aforementioned problems related to large and potentially 
incalculable GUT-scale corrections are automatically avoided in the models we considered 
in Sects. 13.31 and 13.41 All three gauge couplings remain within the perturbative domain 
throughout the entire energy range from irtz to Mqut- 

5 Conclusions 

In this paper we have shown that models with extra Yukawa couplings have a dramatic 
impact on the prediction for a^ at the weak scale. They can bring the two-loop predic- 
tion in line with experimental data without appealing to large GUT-scale or weak-scale 
threshold corrections. This is an effect which has no analogue in the realm of MSSM 
physics - even the top Yukawa coupling is negligible in this context. 

We introduced our main ideas using a simple model with strong GUT coupling. This 
model contains an extra 10+ 10 pair with top-like Yukawa coupling to the MSSM Higgs 
doublets, together with further 5 + 5 pairs making the GUT-scale gauge coupling strong. 
In this context, the main features and implications of our construction can be understood 
in a completely analytical manner. 

We then demonstrated that our promising initial results retain their validity in sit- 
uations where the gauge couplings do not leave the perturbative domain. This more 
complete and partially numerical analysis revealed, among other things, that large input 
values for the extra Yukawa couplings at the GUT scale can lead to an almost perfect 
prediction for the electroweak-scale strong coupling. 

We also tested the sensitivity of our results to extensions of the minimal setting 
described above. In particular, the positive effect of the extra Yukawas on the unifica- 
tion prediction for 0*3 is significantly further enhanced for a complete vector-like extra 
generation with both up-type and down-type Yukawa couplings. In such models, perfect 
agreement with experiment is achieved without invoking excessively large values for the 
extra Yukawa or gauge couplings at the GUT scale. By contrast, introducing several 
copies of 10 + 10 with corresponding Yukawa couplings leads to a weaker enhancement 
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of the Yukawa effect. This can be traced to the increased scale M at which the extra 
multiplets decouple. Such a raised decoupling scale, which is necessary to keep the GUT 
coupling moderate, partially compensates for the positive effect of further extra Yukawas. 

The presence of additional vector-like matter (at comparatively low-energies) and of 
extra Yukawa couplings, which have large values in the ultraviolet, can clearly have a 
significant impact on the MSSM phenomenology. We expect that the enhanced Higgs Z 
factors, which are at the heart of the effect we analyse, would in particular affect the 
values of the Higgs mass parameters m 2 H , m 2 HD , fi 2 and B/j, which are expected in any 
concrete SUSY breaking model. We note that effects on the Higgs mass bounds in related 
non-supersymmetric models have recently been analysed in |46j. 

Finally, as we have discussed in some detail in our paper, the increased value of the 
GUT gauge coupling does not lead to a precision loss of the unification prediction for a 3 . 
The basic reason is that holomorphicity forbids the dominant higher-loop effects, which 
arise only via the non-holomorphic Z factors of MSSM chiral multiplets. However, non- 
perturbative GUT scale corrections clearly need to be controlled. This is easily possible 
in our setting, since it is not necessary to move to the actual strong-coupling point. 
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Appendix A: Gauge group Z-factors in the MSSM 

From the general formula ([TBI we obtain: 

2C F,1 _J_ 

— = — = 0.96 => \nZ u{1) = -0.04 

yi{m z )J \o>i{m Z )J 



Zu{\) 



2C F,2 _3 

2 



«GUT \ b ' 2 I «GUT 

a2(m z )J \a 2 (m z 



Zsum = — P^ = ^P^t = 0.75 =► \nZ su{2) = -0.29 (49) 



«GUT \ b3 _ I QGUT 

a 3 {m z )) \a 3 {m z ) 



Zsu(3) = ;\ = 7 \ = 0.39 =► \nZ sum = -0.93 
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Here we have used a one-loop GUT coupling 27t/q;gut = 27r ■ 24.3 = 153. The quadratic 
Casimirs Cp,i as well as the one-loop MSSM /9-function coefficients are given by: 

w Fti = (^4 1) and h = (f ' i; - 3 ) • (5o) 
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Appendix B: Effect of top in the MSSM 

In this appendix we analyze the effect of the top Yukawa coupling on the two-loop MSSM 
a3-prediction. We recall that in the presence of Yukawa couplings the matter field Z- 
factors factorize as Z = Z G x Z Y , where, as before, the gauge part is given by Eqs. (12) 
and (13). In order to obtain the Yukawa sector contribution we start with the one-loop 
RGEs 

d\nZ Y dlnZfr d\nZ Y , s 

2vr — gg = -3a t , 2tt — -^ = -2a t , 2vr Qi = -a t , (51) 

at at at 

where we have defined a t = y^/Air, with y t being the canonical top Yukawa coupling. 

The subscript '3' in 'C/3' and 'Q3' refers to the third generation quark-antiquark particles. 

We also have the following one-loop equation for a t : 

1 dlna * a 16 13 

27r — - — = boi t — — a 3 — 6a 2 «i • (52) 

at 3 15 

Following [6], we combine Eqs. (15 ip and ( 1521) . finding 

rflnZ^ 8 3 13 1 dlna t 

27r — - = a 3 a 2 oli — 2tt — . (53) 

dt 3 2 30 2 dt y J 

If we now express each of the three gauge couplings through their respective one-loop 

RGEs, 

dinar 1 h 

(54) 



dt 2tt ' 

we can integrate Eq. (153]) analytically: 

13 _3_ _8_ 

7 y , / «GUT V ' 1 / "GUT \ 2b2 ( «GUT \ 3i>3 / «t(M GUT ) 



Z]. ( mz ) = -^Hi_ ^i^ ^^ ^ rv " u y . (55) 

\ ai (mz)/ \^2{mz)J \a 3 (m z ) J V OLt[m z ) J 

As already mentioned, we focus on the moderately large tan /3 region, where m t = ytV 
with v = 174 GeV. Using the low-scale values m t = 173 GeV and yt(mz) = 0.99, we solve 
Eq. (152]) numerically, using the explicit one- loop formulae for aii(t). The resulting high- 
scale value ^(Mqut) = 0.57 is then used to obtain lnZjy (mz) = 0.74. Furthermore, 
Eqs. (ED imply Z Y 3 = (Z Y Hjj f^ and Z Y Qz = (Z^) 1 / 3 . Thus, we find 

Ato P g) = § In Z Y Qz - H In Z£ + A In Z^ = \ In Zfc = 0.32 . (56) 

While, as is well known, this helps in lowering the a 3 prediction, the effect is far too 
small. 



Appendix C: Matter field Z-factors for n>l 

The numerical values of the wavefunction renormalization factors associated with the 
three gauge groups are listed in Table 

Using these results it is then easy to calculate the two-loop corrections to 27r/a 3 
arising from the different matter field Z-factors (cf. Table [7]). 
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Table 6: Gauge group Z-factors 



n 


In Z S u(3) 


In Zsu(2) 


In Z u(1) 


1 


- 1.06 


- 0.33 


- 0.04 


2 


- 1.27 


- 0.38 


- 0.05 


3 


- 1.63 


- 0.46 


- 0.06 


4 


-2.53 


- 0.66 


- 0.08 


4.45 


-3.87 


-0.93 


-0.11 



Table 7: Matter field contributions 



n 



27 
11 



In Zr 



|lnZ, 



-flnZ E 



-ffln^ 



-Sln^ 



11 



In Z H[7 + 



^Z Hd 



1 


-0.72 


-6.28 


0.21 


3.64 


2.07 


-0.46 


2 


-0.83 


-7.46 


0.26 


4.37 


2.49 


-0.54 


3 


-1.01 


-9.44 


0.27 


5.64 


3.19 


-0.66 


4 


-1.42 


-14.40 


0.41 


8.59 


4.95 


-0.94 


.45 


-2.00 


-21.67 


0.57 


13.07 


7.56 


-1.34 
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